Effects of long-range links on metastable states in a dynamic interaction network 
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We introduce a model for random-walking nodes on a periodic lattice, where the dynamic in- 
teraction network is defined from local interactions and E randomly-added long-range links. With 
periodic states for nodes and an interaction rule of repeated averaging, we numerically find two types 
of metastable states at low- and high-£' limits, respectively, along with consensus states. If we apply 
this model to opinion dynamics, metastable states can be interpreted as sustainable diversities in 
our societies, and our result then implies that, while diversities decrease and eventually disappear 
with more long-range connections, another type of states of diversities can appear when networks 
are almost fully-connected. 
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Complex systems usually consist of interconnected 
parts, forming complex networks, and typically exhibit 
complex and unpredictable behavior. While the con- 
cept of complex networks can be applied to various re- 
search areas, physicists have been trying to find gen- 
eral frameworks for the network theory jTH3], and their 
effort can be summed up in three categories. First, gen- 
erating mechanisms of networks satisfying certain sta- 
tistical properties as in small-world (SW) and scale-free 
networks [IHS], and statistical measures on their struc- 
tures have been studied. Second, if nodes have states, 
dynamic behavior of nodes on given networks has been 
studied: for example, opinion dynamics [THTO] and syn- 
chronization problems [TTJ [12] on SW networks. Finally, 
dynamic behavior of both nodes and links for "coevolv- 
ing" networks [I3HI5] . where nodes and links evolve, in- 
fluencing each other, has been studied, especially when 
time scales for node and link dynamics are comparable. 
Here we present a simple model of coevolving networks 
that can mimic mobile agents in finite spaces [T^ ITOl [T7] . 
With additional long-range links, it produces two types 
of metastable states[ni[Tni[in], and we focus on their dy- 
namic behavior by finding their mean lifetimes. 

For our model, we start with a network defined in Ref. 
|T8]with minor modifications. It considers N nodes resid- 
ing on a two-dimensional (2D) periodic lattice {X x Y) 
(see Fig. [TJ. Each node i (1 < i < N) has its location 
{xi, Hi), where < Xi < X and < Hi < Y . At each time 
step t (= 0,1,2,...), all nodes move independently us- 
ing 2D random walks, and interact only with neighbors. 
We call this network as a dynamic interaction network, 
and it is determined by (i) how neighbors are defined, 
and (ii) how nodes move on the lattice. Here we use the 
Moore neighborhood; hence a node is a neighbor of node 
i if that node resides in any of the surrounding 8 loca- 
tions or the current location of i (see the shaded area in 
Fig. [ij . And nodes move randomly and independently in 
four directions (east, west, south, and north) with equal 




FIG. 1. (color online). An illustration of a 2D periodic 
lattice. Neighbors of node i are in the shaded area, and blue 
lines represent dynamic local interactions, while green lines 
represent randomly-added long-range links. 

probabilities. This network is dynamic, and we will re- 
strict ourselves to the case with X = 100 and Y = 10, 
and N = 10^ (= XY), guaranteeing a giant component 
of connected nodes. 

To this network, we add long-range bidirectional links 
(or shortcuts). At t — 0, E time-invariant links are added 
randomly, meaning that every possible pair of nodes has 
the equal probability of getting an additional link. The 
number of all possible pairs of nodes is N{N — l)/2 
(= E„iax), and we use a normalized value e — E/Emax 
(0 < e < 1) as a control parameter. The network has 
two types of links: one representing short-range dynamic 
interactions, and the other representing randomly-chosen 
long-range static interactions. When both types of links 
exist for any pair of nodes, they are regarded as one link. 
We can use the N x N adjacency matrix K(t) to repre- 
sent this network, where Kij{t) is 1 when there exists a 
link between i and j at time t, and otherwise. 

When e = 0, the network has only local interactions, 
and, when e = 1, the network is fully-connected. We can 
look at some of the basic network properties of this net- 
work: the average degree k(e), the clustering coefficient 
C(e), and the average shortest path length L{e), when all 
nodes are uniformly distributed on the lattice[19]. We 
numerically find these values by averaging over 10^ re- 
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FIG. 2. (color online). For A'^ uniformly distributed 
nodes when {N,X,Y) is (1000,100,10), (4000,200,20), and 
(9000,300,30), we observe (a) C(e) and (b) L(e) versus e. In 
(a), we divide e into three regions for the case oi N = 10^. 



alizations: for e = 0, fc(0) ~ 8.99, C(0) ~ 0.604, and 
L(0) ~ 25.8. The probability of having more than one 
component is noticeable when e is close to 0, and only 
the biggest component was used when calculating C(e) 
and L(e) in those cases. In Fig. |2j we show C(e) and 
L(e) [fc(e) ~ 9 + lO'^e for e 1], and we see some re- 
semblance with the classic models of SW networks [H [S]. 
Here the set of e can be divided into three regions: (I) 
C(e) does not change much while L(e) decreases quickly 
(dominated by short-range links with the SW phenom- 
ena); (II) C(e) starts to decrease and reaches the smallest 
value; (III) C(e) starts to increase before reaching 1 at 
e = 1 (dominated by long-range links). Meanwhile, L{e) 
monotonously decreases to 1. The overall behavior is the 
same for cases with different A^'s. 

We then introduce a periodic state (f>i (0 < 4>i < 1) 
to each node i, and observe dynamics on the network de- 
fined above (0^ can be an opinion of an agent i in opinion 
dynamics [IH], or 2Tr(f)i can be a phase of an oscillator i in 
synchronization problems [111 IT^ [201 121) ). For changes of 
0i's, we use the rule of repeated averaging as below. At 
time t, all nodes update their values in parallel using 

c^^it + 1)^ Mt) + ^7^^ E K^jit)A,,it) (mod 1), 

j 

(1) 

where tr is a coupling constant (0 < cr < 1), ki{t) is 
the degree of node i, and Aji{t) represents the difference 
between (f>j and 4>i at time t. We introduce a function to 
represent the difference between two (p values, 

)j -(b., + 1 (if cbj ~cb^< -0.5), 



A(<^j,0, 



< 

/.J -0,-1 (if >0.5), (2) 

t)j — (pi (otherwise). 



and = A(0j , (/j^). Note that Eq. ([T| is the same as the 
one given in Ref. 1 181 with different notations, except that 
there is no threshold for A^i, and that, if we exchange this 
A-function with the sine function, it basically becomes 
the Kuramoto model [22]. used in many synchronization 
problems. 

In a special case with e = 0, it was already found 
that periodic metastable states (along the x direction) 



can spontaneously emerge from random initial conditions 
especially when X <C F, due to periodicities of both the 
state variable and the lattice structure [18]. We call them 
"period-n" states, and they have finite mean lifetimes 
Tn ■ When e ^ 1, t„ for low n is too long to ob- 
tain numerically, but we will see that, as e increases, 
T,!^^ decreases because long-range links make periodic 
metastable states decay faster and eventually disappear. 

As e increases further, however, another type of 
metastable states starts to appear, and they are group- 
ing metastable states, where nodes are divided into m 
separate groups of different (/)'s, independent of their po- 
sitions on the lattice. For these states, numbers of nodes 
in groups and ^-distances between groups should satisfy 
certain conditions: for example, three groups with sim- 
ilar sizes at (/) ~ 0,1/3, and 2/3, respectively, or ~500 
nodes at ~ 0.5, and two other groups of ~250 nodes at 
(f> ~ 0.1 and 0.9, respectively. We call them "m-group" 
states, while m is an odd number |23|. Their mean life- 
times increase as e increases, and as e approaches 1, m- 
group states with higher m's become more stable. Mean- 
while, converged states, where state variables of all nodes 
converge to one value, are absorbing states in every case 
{consensus in opinion dynamics, or synchronized states 
in synchronization problems) . If either the state variable 
or the lattice is not periodic, there is no metastable state, 
and the system reaches a converged state with repeated- 
averaging mechanisms, while e and a only determine the 
convergence time. Converged states can be also called as 
either period-0 or 1-group states. 

Here we numerically study period-n states in detail, 
as well as 3-group states. Usually, period-n states decay 
into period-(n — 1) states, and later into period-(n — 2) 
states, until they eventually reach converged states — this 
is a decay chain. To investigate this decay process closely, 
one needs a way to find what period a given state is in. 
We introduce v for this purpose. First, we divide the 
lattice into equal-sized regions along the x direction. 
(We assume X is a multiple of iV^, and use — 50 
here.) Second, for each region r (0 < r < AV), we find 
the average of state variables, (pr, of nodes inside that 
region. Since the state variable is periodic, care has to be 
taken; for example, the average of 0.1 and 0.9 should be 
instead of 0.5 [A-function in Eq. ^ should be used]. It is 
only meaningful when the range of all values in a region 
is less than a certain value like 0.5, which is true in most 
cases when e ^ 1, because of the repeated averaging 
between neighbors. Once we find all (prS for all r, we 
can finally define v as 



^ A(0r+l,0r) 



r=0 



(3) 



where ^^r^ = 0o- For period-n states, v is close to n, and 
we can observe behavior of a state by observing the time 
evolution of v. 
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FIG. 3. (color online). Time evolution of u for a state as e 
increases from 0.0002 to 0.002 when a — 0.1. Bottom: Scatter 
plots for e = 0.001 at t = (random initial condition), 10 000 
(period 2), 20 000 (period 1), and 30 000 (period 0). 



In Fig. [3] we observe how a random state evolves for 
different e's when a = 0.1. When e < 0.001, the state 
quickly becomes a period-2 state, and it does not decay 
before i = 30 000. But when e = 0.001, we can observe 
the period-2 state decays into a period-1 state, and then 
a period-0 state. Note that lifetimes found here are only 
for this run, and mean lifetimes for e = 0.001 are shown 
in Fig. 5. As e increases further, period-n states tend 
to decay faster, and for higher e's, period-n states are 
reached only for a short time or bypassed (when e = 
0.002, the state tries to reach a period-1 state, but fails 
before becoming a period-0 state). 

When cr is not too small, the state quickly reaches any 
period-n state or a converged state given an initial con- 
dition, and distributions of these states depend on all 
parameters |18j. When mean lifetimes are short enough 
to be measured, these distributions visibly change as time 
evolves, and, if we use the term Pn{t) to represent the ra- 
tios of period-n states to all runs, we can predict that all 
period-n states decay and po(t) eventually converges to 1 
in a reasonable time. In Fig.Bl we use 10'^ random initial 
conditions (one run each), and observe po{t) as t evolves 
when (7 = 1 and 0.1. We see similar behavior for both a 
values. Even though decays from n to n — 1 are all ex- 
ponential, curves don't have to be straight here, because 
period-1 and period-2 states coexist initially and period- 
2 states go through period-1 states to get to period-0 
states. 

If we specifically choose an initial condition as a period- 
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FIG. 4. (color online). Observing po{t), the ratio of the num- 
ber of converged states at t to 10'' runs (each with different 
random initial condition), for several e values when (a) cr = 1, 
and (b) a = 0.1. 
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FIG. 5. (color online). Behavior of tI^^ (1 < 1 < 5) and 
as e varies, (a) rjf^ when cr = 1, (b) when a — 0.1. 
For each parameter set, we used three different realizations 
of period-n states with 10^ runs each. Bars represent mini- 
mum and maximum values. (Insets: From 10^ random initial 
conditions, grouping metastable states appear spontaneously 
as e gets close to 1. In both cases, po is less than 1 when 
e > 0.1.) (c) r{''' and r]®' for < e < 1 when cr = and 0.1. 
There are three regions of e, and the middle region (shaded 
for the case of cr = 1), in which metastable states cannot be 
sustained, shrinks as cr decreases. 



n state and observe the behavior of with 10^ 

runs, then curves are straight, and we can find Tn^'' by 
finding slopes using the method of linear least squares. 
In Fig. g we find r^^^ (1 < n < 5) and the mean 
lifetime of states with three equal-size groups, as e is 
varied when a — 1 and 0.1. We found several interesting 

properties for rlf^ in Figs. 51 a) and 51 b). For period- 

■ ■ • (p) 

n states for low n, there exists a critical e value, en , 

after which period-n states cannot be sustained — for ex- 
ample, e^^^ ~ 0.002 when a = 0.1 (see Vig.^. When 
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e > e„ , period-n states cannot be sustained for even 
a short period of time, while the variabihty of t^'^ with 
different reahzations becomes greater. When e < e„ , 

(n) 

Tn does not depend much on initial conditions as long 
as they are period-n states. Higher-period states can 
be more stable than lower-period ones for certain e's 
(ri^) > 4''^ > ri^) > T^"^ > Tg^P) for a = 0.1 and 
e = 0.001), but we can predict that the order should 
become normal as e approaches 0. On the other hand, 
when e gets close to 1, grouping metastable states can ex- 
ist. Insets prove their existence: when e > 0.1, po from 
10^ runs (each with a different random initial condition) 
can stay less than 1 for a considerable amount of time. 
The closer e is to 1, the longer the mean lifetimes of these 
metastablestates become. 

In Fig. pfc), we observe t[^^ and Tg^-* for the whole 
range of e when a = and 0.1. There also exists a criti- 
cal value Cg^^ for 3-group states. When e > Cg^^ 3-group 
states can exist and Tg®-* increases with e. When e < eg^', 
3-group states cannot be sustained, but Tg^'' tend to co- 
incide with t[^\ increasing again as e decreases, because 
3-group states try to form period-1 state before converg- 
ing to 1-group states without success. Based on results 
for a — 1, we divide e into three regions, which match 
well with those found in Fig.[2|a): (I) < e < e^^'', where 
periodic metastable states can exist (dominated by short- 
range links); (II) e^f'' < e < e'f\ where only converged 
states are sustained; (III) e'f^ < e < 1, where group- 
ing metastable states can exist (dominated by long-range 
links). It is also interesting that the range of e, where 
small-world phenomena are observed, coincides with the 
range where periodic metastable states have measurable 
mean lifetimes. We also observe results for a = 0.1, 
and they show similar behavior with longer lifetimes and 
smaller convergence regions. If a is much smaller, the 
system can behave in an unexpected manner; for exam- 
ple, even when e = 0, grouping metastable states can 
appear (not shown here; to be discussed in future work). 

In summary, we studied dynamics of random-walking 
nodes with periodic states in a periodic lattice, varying 
the number of long-range links, e, and observed two types 
of metastable states in the system at low- and high-e lim- 
its, respectively. Metastable states in our model can be 
interpreted as dynamic, yet sustainable, diversities abun- 
dant in our societies [23], and our result implies that ever- 
increasing connections of our modern societies can make 
existing diversities disappear, but if we are connected too 
well, another type of states of diversities can appear. 

Metastable states can exist due to periodicities of both 
the node state and the space, and we can modify or ex- 
tend this simple model without changing the overall be- 
havior to suit realistic needs. One coupling constant a is 
given to both short- and long-range interactions here, but 
we can have two coupling const ants [12], or it can be het- 



erogeneous for all nodes. We can add noise in Eq. 0, be- 
cause noise can play an important role in clustering[24j. 
In addition, long-range links can be rewired based on 
node states, making the model truly coevolving. 
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